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In models where lepton number is considered to be a spontaneously-broken global symmetry 
a massless Goldstone boson, the majoron (J), appears. We describe a procedure to explore the 
muon-electron-majoron coupling using the results from n — e conversion search experiments. To 
accomplish that, we determine how the energy spectrum of the muon decay into an electron and a 
majoron is modified by binding effects in a muonic atom. We find that the future (x e conversion 
experiments may be able to produce bounds on the n — > ej rate which are comparable with the 
present ones from direct searches. 



I. INTRODUCTION 



The observation of neutrino oscillations has established that at least some of the neutrinos have mass [l[ . The origin 
of the small neutrino mass differences and pattern of mixing angles which are seen in the experiments is not known at 
present. Exploring the nature of neutrino masses and mixing angles may allow us to glimpse particles and interactions 
beyond the Standard Model Q. One class of models that generates neutrino masses, which could provide a possible 
' explanation of the origin of the observed oscillations, considers lepton number as a spontaneously-broken global 
^p" 1 , symmetry. In that case a massless Goldstone boson, the majoron (J) appears. Models of this kind were explored well 
O before neutrino oscillations were experimentally established Q . Some of those models 0, Q predicted large additional 
contributions to the invisible decay width of the Z boson and were excluded after the Z-width measurement at LEP. 
It is possible, though, to construct supersymmetric models with spontaneous breaking of i?-parity (and thus of lepton 
number) where no significant invisible Z-width isjpresent, and the LEP bounds can be evaded |6j. 
£\J • In this paper, we focus on the model of Ref. Jo] which permits charged-lepton decays with majoron emission 

Those processes were recently revisited in Ref. [8J, where it was shown that the /i — > e J decay rate is allowed to be 
large, and could potentially be in a region where it could be measured. It is quite interesting to study those decays, 
since they can explore regions of the supersymmetric parameter space that are not probed by collider searches. 
Experimentally it is quite a difficult task to improve on the current limits for the branching ratio B(fi — » eJ). In this 
• ■ paper we show that the future fi — e conversion experiments [9] may be able to produce bounds that are comparable 
to the present ones. Previous studies of (/i _ ,e + ) conversion mechanisms involved majorons [Toj in now disfavored 
models where the majoron is a gauge non-singlet. 

The current limit for the branching ratio of a muon decaying to an electron and a majoron is [ll| 



B{n eJ) = ^ eJ) _ < 8.4 x 1(T 6 . (1) 

We will study fi — e conversion experiments, which can reach high sensitivities, and possibly improve on this limit. 
The fi — e conversion experiments produce muonic atoms (by stopping muons in a target) and then search for the 
process 

(j,- + (A,Z) -> e" + (A,Z), (2) 

where (A, Z) represents a nucleus of atomic number Z and mass number A. The signal in conversion experiments is 
a mono-energetic electron at energy E^ = — Eb — E lec (with m M the muon mass, Eh the binding energy of the 
muonic atom, and E rec the nuclear- recoil energy). If we now consider the process of /i — > e J decay in the orbit of a 
nucleus (A, Z), i.e. 

lf + (A,Z)->e- + J+(A,Z), (3) 

the outgoing electron can have energies up to the value of the conversion energy E^ e since the nucleus can absorb 
momentum. E^ e is also the maximum electron energy for muon decay in orbit (DIO), which constitutes the main 
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physics background source in the search for the conversion process in Eq. ([2j . Since the majoron J is not observed, the 
signal for the process in Eq. ^ would be electrons below E^ e . Therefore, the measurement of the electron spectrum 
close to E^e, as is done in the conversion experiments, can be used to obtain a limit for B(/i — > eJ). The electron 
spectrum for the decay in Eq. is a delta- function like shape around E e ~ m M /2 (where E e is the electron energy) 
with a tail due to bound-state effects. The spectrum for the two-body free decay [i — > e J is just a delta function at 
E e ~ rn^/2. Even though we only look for events at the high-energy tail of the electron-energy spectrum in /i — > e J 
decay in orbit, the fact that the current limit in Eq. (JT|) is not very stringent, combined with the high sensitivity that 
the future conversion experiments are expected to reach, might allow them to improve on the present constraints. 
The \i — > e conversion experiments provide a limit on -B(/i — > e J) by first determining 

T(^A,Z)^e(A : Z)) = R ^ 

-^capture 

where r capt urc is the rate of the nuclear muon capture, r captU rc = T (/z + (A, Z) — > (A, Z — 1) + v^). This result places 
a bound on the majoron emission in muon decay according to 

Tjn -» eJ) x fj 

N R R,j. e , (5) 

1 capture 

where fj is the fraction of fi — > e J decay in orbit events in the signal region of the conversion experiment and N R is 
a correction factor for the phase space region used in the search for p. + (A, Z) — > e + J + (A, Z). The limit on the 
branching ratio is then given by 

B{» -> eJ) = ^ eJ) _ ~ N f^ T ^-_ . (6) 

In the near future, the DeeMe Collaboration [l2j has proposed to reach i? Me ~ 10~ 14 sensitivity. Furthermore, the 
planned conversion experiments, Mu2e at Fermilab [lJl and COMET at J-PARC aim for sensitivities at the 
R^ e ~ 10~ 16 level; in addition, both experiments are speculating on a later phase which would aim for a sensitivity 
of 10 -18 . Comparing Eq. (|6]) with the present limit in Eq. ([1]) we can examine the potential of those experiments to 
improve on the present bounds, as long as the factor (fjT(fi — > eiy^u e )) / (T ca pt UTe N^) is not much smaller than 1CP 11 . 

The paper is organized as follows: in Sec. [TT] we compute the electron spectrum for the fi — > e J decay in atomic 
orbit. Then in Sec. IIII1 we present our numerical results and discuss in detail the bounds that we can obtain for the 
branching ratio B(pi — > eJ). We conclude in Sec. IIV1 



II. MAJORON EMISSION IN ORBIT 



Our aim is to calculate the electron spectrum for majoron emission /i — > eJ in orbit (MEIO), i.e. the process in 
Eq. ©. 

We write the interaction that mediates majoron emission in muon decay \x — > e J as 

C = figiP R eJ + jxgiP L eJ, (7) 

where Pr = (1 + 7s)/2 and Pl = (1 — 7s)/2. The couplings g\ and gi are dimensionless, i.e. they can be written in 
terms of parameters of the underlying supersymmetric model [§[. The matrix element for the free [i — > eJ decay is 
given by 

M=u(p e )(g 1 P R + g 2 P L )u(p IM ) ; \M\ 2 = u(p e ) {g x P R + g 2 P L ) u( Ptl )u( Ptl ) ( gi P L + g 2 P R ) u(p e ), (8) 

where and p e are the momentum of the muon and the electron, respectively. Summing over the electron spin and 
averaging over the muon spin we obtain 

Y)M\ 2 = {gl+g 2 2 )p e -p^ (9) 

where we considered the electron to be massless. We always consider a massless electron in this paper, since we are 
interested in the high-energy part of the spectrum for MEIO. The free decay rate, To, is then given by 

r » = *b / A p|k (2 - ,45<4,(p " - - - w »s^i 2 
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9l + 9l f El / _ m, \ {gl + g\)m 



167T ./ TO^ — i? R 



.j*^ (10) 



where p,j — (Ej,pj) is the majoron 4- momentum. Obviously, since this is a 2-body decay, the electron-energy 
spectrum is just a 5 function at E e — For MEIO Eq. (fTU|) gets replaced by 

e spin 

where = m M — i?t>, and 

J:= f d\e-^ ^ e { 9l P R + g 2 P L )tPn, (12) 



where ip e and <^ are the solutions of the Dirac equation (in the potential created by the nucleus) for the electron 
and the muon, respectively. We incorporate the average over the muon spin in the definition of ip^, while we do 
not incorporate the sum over the electron spin in the definition of (p e . When the muonic atom is formed the muon 
cascades down almost immediately to the ground state, and this process also depolarizes the muon [l5j]. We consider 
an unpolarized muon in the IS* state and write the muon wavefunction as 



(o=e°.( J' (13) 



where a s is the amplitude of the muon state with spin projection s. For an unpolarized muon we have |a s | 2 = 1/2. 
Xk — Xk(^) are the spin-angular functions, which are given by 



with C(lsj;l z s z j z ) the Clebsch-Gordan coefficients, are the spherical harmonics and x™ are the spin 1/2 eigen- 
functions. G and F are the solutions of the radial Dirac equations, which are taken to be normalized as 

r 2 (F 2 + G 2 )dr = 1. (15) 

The electron wavefunction is expanded in partial waves according to 



Kfl Kfl 



where g K and f K are the solutions of the radial Dirac equations for the electron, labeled by k = ±1,±2, . . . [16(, t is 
the z-component of the electron spin, and the a Kfl t coefficients are given by 

a K , t = i l «-j=C (l K lj K ^-tt^\ Yt- t *{p e )e-^, (17) 

where 8 K is the Coulomb phase shift (the distortion from a plane wave due to the potential of the nucleus), j K — 
\k\ — 1/2, and l K = j K — Sign(«;)/2. The electron wavefunctions are normalized in the energy scale, according to 

d 3 r^ K * w ^,w> = 2n5^5 KK/ S(W - W), (18) 

where ip^ w corresponds to a solution with energy W. 

We ignored nuclear-recoil effects to write Eq. (JTTJ) but will incorporate them later. Integrating over the solid angle 
in Eq. (fT2"j) we obtain 



2 

KM sKfj, 



x {(gi + .92) (g K G - f K F) 5 KIk + i( 9l - g 2 ) (f K G + g K F) S Kl _ K } . (19) 
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We can then write the result for the spectrum, 

1 dT 



T Q dE e 



E 

Kk 



irm, 



-(2j K + l)[E fl -E e ]\S KK \ 



where 



Skk 



drr 2 j K QEp - E e ]r) {{g K G - f K F) 5 K i K + % (f K G + g K F) 6 m _ K } 



(20) 



(21) 



In Eq. (f2T)[) . the sum over K goes from to oo, and for a given value of K, k can only take the values ±K and 
±(K + 1), but k can never be equal to 0. j n {z) is the spherical Bessel function of order n. 

Since future \i — > e conversion experiments may use aluminum or heavier nuclei as targets, the nucleus is, at least, 
more than 200 times heavier than the muon and nuclear-recoil effects are negligible for most of the electron spectrum. 
Recoil effects could be important close to the high-energy endpoint since the y m odify the maximum allowed electron 
energy. In that region we can approximate the nuclear-recoil energy E rec as [171 ] 



Ever 



E, 



\VN?_ 

2mjv 2m 



N 



(22) 



with tun the mass of the nucleus and pn its three- momentum. Within this approximation, the inclusion of recoil 
effects reduces to the substitution 



Efjt — E e 



Efj, — E e 



11 
2m 



N 



(23) 



E u — E 2 J (2m n). Inclusion of nuclear- 



inside the square brackets in Eqs. (|2"0|) and (l2~Tj) . The endpoint is now at E e — ^ , 
recoil effects beyond the approximation in the equations above is unnecessary for all practical purposes. Our final 
result for the electron spectrum in MEIO then reads 



1 dT 
Tn~dE~, 



E 



i 



-(2j« + l) 



Kk 



E,, — E P 



El 
2m n 
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Kk 



with 



Skk '■= j 



drr j K 



E n — E P 



El 
2m n 



-) {{9kG - f K F) 5 KIk + i {f K G + g K F) 5 m _ K } 



(24) 



(25) 



A. Taylor expansion around the endpoint 

We can perform a Taylor expansion of the electron energy spectrum in Eq. (|24p around the endpoint to make its 
behavior manifest. For that, we need the Taylor expansion of Sk k m Eq. ([23)) . We obtain, noting that the electron 
wavefunctions g K and /« depend on E e , 

Skk = J drr 2 (g_iG-/_iF) S 0K 5-i K 

+ f drr 2 ^- + [i(f-iG + g-iF)S 1K 5- lK + (g-2G- f- 2 F) S 1k S- 2k ] 

- (g'-iG - fL x F) 5 0K 5^ lK } (e„ - E, 



where 



ri _ df R , _ dg K 

dEp ; 9k '- dEp- ( <) 



It is understood that the electron wavefunctions in Eq. (|26[) (or their derivatives) are evaluated at the endpoint energy 
E e — E^ — E 2 /(2mpf). The corresponding Taylor expansion without including recoil effects is recovered by taking 
the limit m^r — > oo. Since we have the property 

drr 2 g_ 1 G= f drr 2 f_ ± F, (28) 
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TABLE I. Values for the parameters 


in the Fermi distribution in 


Eq. ([291) 0. nuclear 


masses, muon enerj 


U Efj,, and endpoint 


energy _B Me , 


for the elements used in 


the text. 








Nucleus 


ro (fm) 


a (fm) 


m N (MeV) 


£ M (MeV) 


(MeV) 


Al(Z = 13) 


2.84 


0.569 


25133 


105.194 


104.973 


Ti(Z = 22) 


3.84 


0.588 


44588 


104.394 


104.272 


Au(Z = 79) 


6.38 


0.535 


183473 


95.533 


95.508 



when the energies of the muon and the electron are equal 17], the leading term in the expansion of Sk k vanishes 
when we neglect nuclear-recoil effects. Therefore, the Taylor expansion of the electron spectrum without nuclear-recoil 
effects starts at order (E^ — E e ) 3 . If we take nuclear recoil into account the leading term of the Taylor expansion for 
Skk no longer vanishes exactly, but it is suppressed by the inverse of rajy- 



III. BOUNDS ON THE BRANCHING RATIO 



Equation (IM1) gives the electron spectrum for MEIO. We will use it here to obtain bounds for the branching ratio 
B(p — > eJ). We solve the Dirac equation numerically to obtain the electron (g K and /«) and muon (G and F) 
wavefunctions that appear in Eq. (I24[) . To do that, we consider a nucleus of finite size, characterized by a two- 
parameter Fermi distribution p(r), given by 

p(r)=Po W- ( 29 ) 

1 + e^ 

The parameters ro and a, the nuclear masses, the muon binding energy E^ 7 and the endpoint energy E^ e for the 
elements of current experimental interest are summarized in Table |TJ The normalization factor po in Eq. (29) is 
determined from the condition J d 3 rp(r) = —Za. For the muon mass and the fine structure constant we use the 
values to m = 105.6584 MeV, a = jg^ggg ■ Radiative corrections were not included in Eq. (|24l) but are not expected to 
significantly modify the results. Note that, as it happens for the usual muon DIO (see Ref. [13); to obtain the correct 
result for the high-energy tail of the spectrum it is necessary to include finite nuclear size effects, the interaction of 
the outgoing electron with the field of the nucleus, the Dirac (rather than Schrodinger) wavefunction for the muon 
and (at least for not very heavy elements) nuclear- recoil effects, as we did. 

Conversion experiments measure the electron spectrum in a window of a few MeV near the energy E^, in search 
for a peak that would reveal the conversion process in Eq. ([2]). In this energy window, electrons from muon DIO arc 
also present and are seen by the experiments. The spectrum for muon DIO has been recently studied in detail [l9[ and 
it is now under good theoretical control. Electrons coming from MEIO would appear as an additional contribution 
on top of the electrons coming from muon DIO. 



A. Existing conversion results 



Currently the most stringent upper limit for the conversion branching ratio is g iven by the SINDRUM II Collabora- 
tion, i? Me < 7 x 10 -13 (90% C.L.), using a gold (Au, Z = 79) target [2(|. Ref. [20j measured the electron spectrum in a 
region from 90 MeV to E^ e . From Eq. (|2~4")l we obtain the electron spectrum for MEIO for gold, which for illustration 
we plot in Fig. [1] as the solid line. We find that the MEIO spectrum for energies E e > 90 MeV is very well fitted by 



where 




Since the gold nucleus is quite heavy, the terms proportional to S and S 2 in the spectrum are found to be negligibly 
small (see the discussion in Sec. Ill A|) . The fraction fj of p, e J decay in orbit events for energies above 90 MeV is 
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FIG. 1. Electron spectrum for majoron emission in orbit for gold (solid line). The second panel is a zoom for E e > 90 MeV. 
The dashed line in the second panel is the electron spectrum for DIO in gold, multiplied by a constant (C = 333) to make it 
coincide with the MEIO rate at E e = 90 MeV. 



given by 



fj\Au,E e >90 MeV = / f-rfP"^ 8 = 2A X 10 ' (32) 



The total muon lifetime in gold is 88 ns, determined by the capture rate. Therefore, using the estimate for the limit 
on the branching ratio that we derived in Eq. ([6]), we have 

n , t \ NftRue r capt ure NrR^ 2.197 fiS NrRub nr /00 n 

B(fi->eJ) 7~ r/ ~f , 55 (33) 

where we also used that the free muon width is T(fi —¥ ev^v e ) — l/(2.197//s), and Nr is the factor rep resenting the 
event distribution observed below the conversion peak. The upper limit R^ e < 7 x 10~ 13 set by Ref. [2(| was obtained 
from a likelihood analysis considering the mono-energetic conversion signal, muon DIO, and additional backgrounds 
from radiative muon capture, pion decays and cosmic rays. The shape of the electron energy distribution seen in the 
experiment was well described by muon DIO. It is reasonable to assume that the events seen in the region E e > 90 
MeV came from muon DIO and that no additional events from majoron emission were present. Nevertheless, to avoid 
doing background subtractions of DIO events, we will consider that the number of MEIO events cannot be larger 
than those seen in the experiment. Six events were seen for E e > 90 MeV, and if we assume a Gaussian distribution 
(with mean and variance equal to 6), the 90% C.L. limit corresponds to 9 events. Therefore, to obtain the limit on 
B(fi eJ), we use Eq. (|3"3"|) with the factor Nr = 9/2.3 (where the 2.3 in the denominator corresponds to the events 
for 90% C.L. on a Poisson distribution with expected number of occurrences 1 or 0). Doing that we obtain 

B{» - e J) < ^ 25 ~ f^ll 3 25 ~ 3 x 10-. (34) 

The upper limit in Eq. (|34p is still less restrictive than the current limit in Eq. (J]). Thus, the existing conversion data 
does not improve on the present limits for the [i — e — J coupling. 



B. New conversion experiments (Mu2e and COMET) 

Future /i — > e conversion experiments Mu2e at Fermilab [l3[ and COMET at J-PARC [l4| aim to reach sensitivities 
at the 10~ 16 level, and to use an aluminum (Al, Z = 13) target. The ratio of muon capture width in Al over the free 
muon width is around 1.5 (as compared to 25 for Au, as in Eq. (|33l) ). but fj is lower due to the much lower Z. The 
fi — > e conversion energy in Al is around 105 MeV (see Table |T]). Using E e > 100 MeV as the signal region, fj for Al 
is 2.2 x 10~ 10 , and Nr = 27 based on subtraction of the DIO background, resulting in 



NrR^ r capture NrR^ 27 x HT 16 

fj T{fi -> e^D e ) ~ fj ' ~2.2x 10- 



B(jjl -> eJ) ~ , T„""V. \ ~ tlJ T^ 1 ^ - o oV/in-io 1 - 5 ~ L9 x 10_5 ' ( 35 ) 



which is comparable to the current limit in Eq. ([T]). If a sensitivity of 10 18 is reached, the limit on -> eJ) will 
improve by an order of magnitude due to improved statistical precision on DIO. 
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FIG. 2. Electron spectrum for majoron emission in orbit for Al (solid line). The second panel is a zoom for E e > 100 MeV. 
The dashed line in the second panel is the electron spectrum for DIO in Al, multiplied by a constant (C = 415) to make it 
coincide with the MEIO rate at E e = 100 MeV. 



To obtain a more detailed estimate of the upper limit for — > eJ) that Mu2e and COMET could obtain, we 
need to convolute the MEIO spectrum with the experimental energy resolution modeled by those collaborations. Once 
again, from Eq. we get the electron spectrum for MEIO for Al, which for illustration we plot in Fig. [5] as the 
solid line. We find that the spectrum for energies E e > 100 MeV is very well fitted bjQ 



i dr 



Al, £ e >100 MeV 



= (3.289 x IQ- 10 5 + 3.137 x 10~ 7 <S 2 + 1.027 x 10~ 4 <5 3 + 1.438 x 10~ 3 (5 4 + 2.419 x 10~ 3 (5 5 + 1.215 x 10" 1 S 6 ) 



m 



(36) 

Convoluting the above equation with the estimated energy resolution of the experiment (which we denote as F(E e )) 
we can obtain a value of fj 



fj\ 



E c >x 



1 dT 

T lE e 



F dE e , 



(37) 



where <g> denotes the convolution. Using the estimated energy resolution for Mu2e and COMET, and estimating the 
expected number of DIO events in the signal region (iVoio — 2400), we find that the bound on B(/i —> eJ) that a 
conversion experiment with R^ e ~ 10 -16 on Al may be able to place will be B{^l — > eJ) < 2 x 10~ 5 , roughly at the 
same level as the current one in Eq. ([1}. 

Since COMET and Mu2e also consider Ti as a viable target we give, for completeness, the polynomial that fits the 
MEIO spectrum in Ti for E e > 99 MeV0 



1 dT 
ro" dEe 



Ti, B e >99 MoV 



= (5.404 x \Q- W 5 + 9.301 x 10~ 7 <5 2 + 5.552 x 10~ 4 <5 3 + 8.113 x 10~ 3 (5 4 + 5.470 x 10~ 2 <S 5 + 4.244 x 10" 1 S 6 ) . 

rrifj, 

(38) 

The bounds on B(fi — > eJ) we would obtain with a Ti target are similar to those for Al. 



C. Discussion 



We have shown in the previous section that a /i — > e conversion experiment with an Al target and sensitivity at 
the 10~ 16 — 10~ 18 level may be able to produce bounds on the [i — e — J coupling which are competitive with the 



1 Eq. (1 36 I t may be used in a wider energy range. It reproduces Eq. 1241 ) for Al with an accuracy better than 10% until E e ~ 90 MeV. 

2 It can be used until E e ~ 93 MeV to reproduce Eq. i ]24fl for Ti with an accuracy better than 10%. 
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present ones. If the electron distribution seen in the experiments around Ea e agrees well with DIO, the procedure 
we have presented allows to place a bound on B(fi — » eJ). If it is found that the electron distribution is not well 
described by DIO, it could be checked if the addition of a component following Eq. (|36|) improves the agreement 
which could be consistent with the presence of events coming from MEIO. Note that (when we neglect nuclear- recoil 
effects) the electron spectrum for MEIO goes as ~ (-E>e — E e ) 3 near the endpoint, while the usual muon DIO goes as 
~ (E^e — E e ) 5 . MEIO is, therefore, less suppressed than DIO in this region, and the endpoint is a favorable region 
to search for the /i — > eJ process. To compare the shapes of the two processes near the endpoint, we also show in 
the second plot in Fig. [5] (Fig. [TJ the electron spectrum for DIO in Al ( Au) as the dashed line [l9[ , multiplied by a 
constant C = 415 (C = 333) to make it coincide with the MEIO value at E e = 100 MeV (E e = 90 McV). 

One may also ask which target materials give better sensitivity to B(fi — » eJ). When we increase Z, the fraction fj 
increases (which improves the bound one would obtain), but the capture width and Nr also increase (which worsens 
the bound one would obtain). The net result is that, for a given sensitivity to the conversion search, the sensitivity 
to B(fi — > eJ) is roughly independent of Z. 

IV. CONCLUSIONS 

We have computed the electron spectrum for the fj, — > e J decay when the muon is orbiting a nucleus. Using those 
results, we described a procedure to probe the muon-electron-majoron coupling using /j — > e conversion experiments. 
Using results of conversion experiments to improve the limit on B([i — > e J) does not require any dedicated experimental 
search. B(pi — > eJ) can also be probed in /i — » ej searches [8|, although some relaxation in the cuts of those 
experiments is required to improve the present limits. The future conversion experiments, Mu2e and COMET, may 
have the capability to produce bounds on B{\i — > e J) that are competitive with the present ones, and possibly improve 
them. The results presented in this paper strengthen the physics case for the upcoming conversion experiments. 
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